Equivalence of critical scaling laws for many-body entanglement in the 

Lipkin-Meshkov-Glick model 
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We establish a relation between several entanglement properties in the Lipkin-Meshkov-Glick 
model, which is a system of mutually interacting spins embedded in a magnetic field. We provide 
analytical proofs that the single-copy entanglement and the global geometric entanglement of the 
ground state close to and at criticality behave as the entanglement entropy. These results are in 
deep contrast to what is found in one- dimensional spin systems where these three entanglement 
measures behave differently. 
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Introduction — In the field of quantum many-body 
physics, common wisdom dictates that quantum entan- 
glement plays a key role in the occurrence of important 
collective phenomena at zero temperature [1]. Under- 
standing the entanglement properties of many-body sys- 
tems at the critical points of quantum phase transitions 
is, therefore, a great theoretical challenge. In this respect, 
considerable efforts have been devoted in recent years to- 
wards a theory of entanglement in extended systems (see 
[2] for a review). 

Within this context, significant attention has been paid 
to systems in one dimension (ID), which are mostly 
tractable by analytical studies. In particular, the pio- 
neering works [3, 4, 5] established that the ground state 
entanglement entropy £ obeys universal scaling laws in 
critical regions (that is, either close to or at criticality) 
that are described by an underlying conformal field the- 
ory [G] . This entropy measures the entanglement between 
two subspaces and thus relies on a bipartition of the orig- 
inal system. Along the same line, universal scaling laws 
have been also obtained for the single-copy entanglement 
S [7, 8] which quantifies the amount of entanglement dis- 
tillable from a single specimen of a quantum system, and 
for the density of global geometric entanglement per sub- 
system Q /N [9, 10] which measures the distance, in the 
Hilbcrt space, to the closest separable state. Surprisingly 
enough, for critical ID systems, these three quantities 
which are a priori very different, turn out to be deeply 
intertwined since they obey Q/N ~ S/2 ~ £/A. This 
result has important consequences in our understanding 
of renormalization group flows in ID [11, 12, 13, 14]. 

It is then natural to wonder whether such a relation 
still holds for higher-dimensional critical systems. How- 
ever, answering this question implies to solve some diffi- 
cult problems. Indeed, the majority of systems beyond 
ID do not admit an analytical solution and their proper- 
tics need to be unveiled by numerical simulations, always 
with a partial success [15, 16, 17, 18]. 

The aim of this letter is to investigate this issue by pro- 



viding an analytical derivation of the relation at critical- 
ity between the entanglement entropy £, the single copy 
entanglement S and the global geometric entanglement Q 
in the Lipkin-Meshkov-Glick (LMG) model [19, 20, 21]. 

Originally introduced in nuclear physics, the LMG 
model has, since then, been used in the description of 
many physical systems among which are two-mode Bosc- 
Einstcin condensates [22] or small ferromagnetic par- 
ticles [23]. However, the full spectrum of this model 
has only been exactly determined recently in the ther- 
modynamical limit [24] and has revealed a very rich 
structure. Some entanglement properties have already 
been investigated in this model (concurrence [25, 26, 27], 
cntropy[28, 29, 30], fidelity [31]) albeit, contrary to ID 
systems, there is no comparison amongst different mea- 
sures. 

In this work, we bridge this gap by computing ex- 
actly, in the thermodynamical limit, the global geomet- 
ric entanglement as well as the single-copy entanglement 
of the ground state. These results allow us to extract 
their behavior in the critical region and to establish that 
these quantities obey exactly the same scaling laws as 
the entropy, i.e., Q ~ S ~ £ in deep contrast with one- 
dimensional spin systems. 

The model — The LMG model describes a system of N 
spins 1/2 mutually interacting and embedded in a trans- 
verse magnetic field. Its dimensionless Hamiltonian is 
given by 



■ 7 S v ) ~ hS z , 



(1) 



where S a = 2&Li a a/^ are the total spin operators in the 
direction a, o~ l a is the Pauli matrix a for spin i, 7 is the 
anisotropy parameter and h is the transverse magnetic 
field. Notice that the Hamiltonian (1) can be seen as 
those of a iV-dimcnsional system of N spins 1/2 or a 
zero-dimensional system of one spin N/2 particle. Here, 
we focus on the ferromagnetic case and, without loss of 
generality, we assume < 7 < 1 and h > 0. 
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It is well-known that the system undergoes a second- 
order quantum phase transition at h = 1 which separates 
a symmetric phase for h > 1 from a broken phase h < 1 . 
The basic properties of these phases can be understood 
in terms of a mean-field approach [27, 32]. Within this 
approximation, the ground state is described by a fully- 
polarized state which is unique for h > 1, and two- fold 
degenerate for h < 1. However, although the mean-field 
treatment perfectly describes the transition, the exact 
ground state is not a product state, even in the thermo- 
dynamic limit, as initially observed for the concurrence 
[25]. This counter-intuitive result can be easily under- 
stood in terms of quantum fluctuations around a classical 
state by using a Holstcin-Primakoff representation of the 
spin operators [27]. 

In the following, for simplicity, we further restrict our 
study to the symmetric phase and the critical point, that 
is, the region h > 1. 

The global geometric entanglement — To introduce 
this measure, let us consider a pure quantum state of 
N parties |*)eW = <g)? =1 , where U [i] is the Hilbert 
space of party i. We wish to quantify the global multipar- 
tite entanglement of |^/). Following [33], this is achieved 
by considering the maximum fidelity A max between Vf) 
and all the possible separable states |$) of the TV parties 

A nuut = max|($|*)|. (2) 

This quantity can be seen as the distance, in the Hilbert 
space, between the state and the closest separable 
state. Here, a state |$) is said to be separable if it can be 
written as $} = <g>^ 1 (z*i| |)i +«t| where | f)i (| 
is the eigenstate of o\ with eigenvalue +1 (—1). Thus, 
we assume that each individual spin constitutes a single 
party by contrast to the analysis done in Refs. [9, 10], 
where the parties are blocks of several spins. In order 
to have a well-defined measure of entanglement (which 
vanishes when |\&) is a product state) one defines the 
global geometric entanglement Q of state \^f) as 

g = -lnAL x . (3) 

Notably, for the LMG model, one already knows the 
state |$) which maximizes the fidelity. Indeed, one knows 
that the ground state is in the maximum spin sector 
S = N/2 which is permutation- invariant. This implies 
that the closest separable state in the Hilbert space must 
be of the form |$) = ®f = i{u\ J \)i + v\ One thus seeks 
for a coherent state which is as close as possible to the ex- 
act ground state. The answer is given by the mean-field 
treatment detailed in Refs. [27, 32] and, in the symmetric 
phase, one gets u = 1 and v = 0, i. e., the fully-polarized 
state in the z-direction. The next step consists in com- 
puting A max which is more involved since one does not 
know the exact ground state analytically. Nevertheless, 
in the thermodynamical limit (large N), as early given 
in the seminal paper [20] , one can obtain a recursion re- 
lation for the coefficients in the Dickc states basis. After 
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FIG. 1: (color online). Behavior of the geometric entangle- 
ment (left panel ) and single-copy entanglement (right panel) 
as a function of the magnetic field for different system sizes 
N = 16, 32, 64, 128 and 256 (from bottom to top). Red lines 
are the thermodynamical limit given in Eqs. (6) and (12). 
Here, Q has been obtained by a numerical minimization of 
the fidelity over all coherent states. 

simple algebra, one gets the following expression of the 
ground state 

N/2 , \ 1/2 / \ i 

l*o) = (l-t 2 ) 1/4 £( 2 /) (f) (4) 

where 

t = 2/i-7- l-2y/(fe-l)(fe- 7 ) (5) 
1-7 

Here, the state \2i) denotes the eigenstate of S 2 and S z 
with eigenvalues y(y + 1) and -j — 2i. The maximum 
fidelity is then directly given by the coefficient on the 
state |0) so that, in the thermodynamical limit, the global 
geometric entanglement reads 

Q{l,h) = - l -\n{l-t 2 ). (6) 

As can be checked in Fig. 1 (left), this is in perfect agree- 
ment with numerics. 

Near the critical point h = 1, this quantity behaves as 

0( 7j h) = -iln(/i-l) + iln(l- 7 )-ln2 

+o[(h-iy/ 2 }. (7) 

Quite importantly, the above relation shows that when 
reaching the critical point Q diverges as — j hx(h— 1), ex- 
actly as the entanglement entropy £ [29] . We emphasize 
this result is completely nontrivial since Q and £ are, by 
construction, very different objects. As we shall now see, 
it is even more surprising to see that the scaling laws 
derived above are also valid for the single-copy entangle- 
ment. 



3 



The single-copy entanglement — Let us now introduce 
the next entanglement measure for a quantum many- 
body state I*) with reduced density matrix p L for a sub- 
set of L constituents. As explained in Refs. [7, 8], the 
so-called single-copy entanglement S between this subset 
and the rest of the system is given by 



InAi, 



(8) 



where Ai is the largest eigenvalue of the reduced density 
matrix p^. 

Our aim is thus to compute and diagonalizc the re- 
duced density matrix obtained by tracing out over 
(N — L) spins in the LMG model. As detailed in Ref. [27], 
the Hamiltonian (1) can be mapped, in the thermody- 
namical limit, onto a quadratic form of a single bosonic 
mode a via a standard Holstcin-Primakoff transformation 
of the spin operators. As early shown by Bombelli et al. 
[34], one also knows that the reduced density matrix for 
cigenstatcs of a quadratic form can always be written as 
Pl = e~ K with 



K = kq + Kia) a + k,2(o) 2 + a?) 



(9) 



where the coefficients kq, k\ and k>2 have to be deter- 
mined by self-consistent relations. Following Refs. [29, 
30], after diagonalization of K, one gets for the LMG 
model 



PL 



pL+l 



(10) 



where g is a bosonic mode that diagonalizcs K, and 

H = a ~ 1/2 y/[ra + (1 - t)][t + a{\ - r)], (11) 



with t = L/N and a = y/Qi - l)/(h - 7). The pseu- 
doenergy reads e = ln^^ty^. The single-copy entan- 
glement of a subset of L spins, in the thermodynamical 
limit, is thus given by 



<5(r,7,/i) = In 



(12) 



and perfectly matches with numerics as can be seen in 
Fig. 1 (right). 

As previously, one can expand S in the vicinity of the 
critical point to get the relation 

S(r,l,h) = -iln(/ l -l) + iln(l- 7 ) + iln[r(l-r)] 

(13) 



In2 + G[(^-1) 1 / 4 ]. 



Remarkably, the leading terms in the above expression 
are identical to those found for the entanglement entropy 
[29], the only difference occuring in the subleading correc- 
tions, which arc O [{h - l) 1 / 4 ] for S and O [{h - 1) 1/2 ] 
for £ (and Q). 
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FIG. 2: (color online). Behavior of the Q and S at the critical 
point as a function of N. The blue dotted line has a slope 
1/6, as predicted by the scaling hypothesis. The inset shows 
the behavior of the coefficient \ (see text) as a funcion of N. 
The dotted line in the inset also has a slope 1/6. 



Finite-size behavior at the critical point — Using the 
scaling hypothesis discussed in Refs. [27, 29], one can 
further obtain the behavior of Q and S at the critical 
point. For the global geometric entanglement, one gets 

g( 7j h = l)~ilnJV+ilii(l-7). (14) 
o 6 

Note that in the present case, the density of global ge- 
ometric entanglement Q/N vanishes in the thermody- 
namic limit contrary to ID systems where it remains 
finite [9, 10]. For the single-copy entanglement, one sim- 
ilarly obtains 

5(T, 7j h= 1) ~ ilnJV+iln(l- 7 ) + iln[r(l-r)]. (15) 
6 6 2 

Here again, in the large N limit, one recovers the same 
finite-size behavior as for the entropy so that we can fi- 
nally formulate the central result of this paper : in the 
critical region of the LMG model, one has 



(16) 



To check the scaling laws (14) and (15), we computed 
numerically the behavior of Q and S at h — 1 and 7 = 
(and r = 1/2 for S) as a function of the system size, in 
the range N = 16 to N = 8192. Results are shown in 
Fig. 2, where the two straight lines have slope 1/6. For Q, 
the asymptotic regime is nearly already reached for the 
maximum size N = 8192 we used, which confirms (14). 
This is not the case for S, because subleading corrections 
are more important. In order to quantify this, we have 
fitted S with the law S = A + x hi N in the vicinity of 
each value of N. The inset in Fig. 2 shows that \ ~ 
1/6 + BN- 1 ^. Thus, x indeed converges towards the 
expected value of 1/6. This convergence is however very 
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slow, which explains why the asymptotic value is not yet 
reached even for sizes as large as N = 8192. 

Comparison with one- dimensional systems — As said 
in the introduction, the above result does not hold in 
ID systems. Indeed, let us call £ the correlation length 
of a ID system and L the size of a given block. Let us 
also define a quantity I such that I = L at criticality and 
I = £ away from criticality and in the regime L>^. In 
critical regions of ID systems that are described by an 
underlying conformal field theory with central charge c, 
one has that [3, 4, 5, 8, 9, 10] 

~ ls(l) ~ \e{l), (17) 

where £(l) ~ (c/3)lni, and N — ► oo is the number of 
blocks of size L that define each party in the case of 
the geometric entanglement [9, 10]. Therefore, in ID 
one has that S(l) ~ (c/6)lnZ and G(l)/N ~ (c/12)ln/. 
Notice the difference between Eq. (17) and our result 
for the LMG model. The appearance of the factors 1/2 
and 1/4 in front of the single-copy entanglement and the 
entanglement entropy in Eq. (17) seems to be endowed 
with the (assumed) conformal structure at the core of 
the critical points in ID. In the case of the LMG model, 
this conformal structure is no longer present and Eq. (16) 
holds instead of Eq. (17). 

Discussion — Our results have a clearcut interpreta- 
tion from the perspective of quantum information. More 
precisely, Eq. (16) establishes that the ground state of the 
LMG model is equally suited as a resource for two dif- 
ferent tasks, namely, the concentration of entanglement 
by local operations [35] from (i) infinitely-many copies 
of the system (quantified by £), and from (ii) just one 
system (quantified by S). As shown in this paper, the 
capability to perform such tasks is also equivalent to the 
fidelity between the ground state of the system and the 
closest separable state of all the spins (quantified by Q). 
Our results also suggest that the relationship between 
entanglement measures depends on the dimension of the 
quantum system considered. An interesting issue would 
be to analyze these measures in exactly solvable systems 
such as, for instance, the celebrated Toric Code model 
[36]. However, note that if the definition of Q may be 
non ambiguous in any dimension, it is not the case for 
other measures. 
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